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Abstract 


Fractal  and  multifractal  methods  are  used  to  examine  the  changing  complexity  of  active  regions  on  the  solar 
disk.  It  is  shown  that  to  accurately  recover  the  fractal  and  multifractal  spectrum  of  the  active  regions,  it  is 
necessary  to  remove  contribution  to  the  spectrum  from  the  surrounding  quiet  Sun  features.  A  segmentation 
method  in  wavelet  modulus  space  is  presented  to  remove  the  quiet  Sun  contributions  from  the  subsequent 
calculations.  The  method  is  shown  to  greatly  increase  the  accuracy  of  the  spectrum  and  allows  for  accurate 
comparison  of  an  active  regions  multifractal  characteristics  and  its  probability  to  flare. 


1  Introduction 
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1  Introduction 


Solar  active  regions  are  highly  dynamic  complex  systems  that  have  been  studied  with  varying  degrees  of 
success  using  the  fractal  and  multifractal  tools.  One  of  the  main  disadvantages  of  current  methods  is  their 
inability  to  distinguish  and  separate  quiet  Sun  features  from  those  of  active  regions.  Quiet  Sun  magnetic 
features  are  shown  to  be  statistically  different  from  that  of  active  regions.  As  such,  methods  which  fail  to 
separate  quiet  Sun  features  from  the  calculation  of  the  active  region  fractal  and  multifractal  parameters  have 
decreased  accuracy.  Using  the  proposed  segmentation  and  multifractal  methods,  a  correlation  is  presented 
between  the  multifractal  parameters  and  active  region  flaring. 

Natural  systems  are  complex  and  chaotic  in  their  behaviour.  As  such,  simple  mathematical  tools  are  unable 
to  describe  them  accurately.  Fractals  were  first  proposed  to  solve  this  issue  in  the  study  of  turbulence  and 
complexity  (Schroeder,  1991;  Mandelbrot,  1977).  Fractals  are  self-similar  objects,  meaning  that  at  each  level 
of  detail  a  similar  pattern  or  structure  is  observed.  Fractals  described  the  geometry  of  an  object  at  each 
of  these  scales  and  are  related  to  an  objects  Euclidean  geometry.  As  such,  lines  and  squares  have  fractal 
dimensions  of  one  and  two  respectively.  An  object  with  a  fractal  dimension  less  than  one  has  less  detail 
than  a  line  and  is  dust-like.  Objects  with  fractal  dimensions  between  one  and  two  have  more  detail  in  their 
structure  than  that  of  a  line  but  do  not  occupy  fully  the  space  of  a  square. 

Multifractal  are  a  natural  extension  to  the  fractal  idea,  developed  to  characterise  systems  with  more  than 
one  fractal  dimension,  as  found  in  chaotic  dynamical  systems  (Mandelbrot,  1974).  Multifractal  objects  have 
a  spectrum  of  fractal  dimensions.  Each  fractal  dimension  contributes  with  varying  degrees  of  importance 
to  the  overall  design  and  structure  of  the  system. 


2  Methods 


Fractals  and  multifractals  measure  the  detail  across  scale  of  self-similar  objects.  In  this  section,  the  tradi¬ 
tional  box-counting  algorithm  is  described  and  associated  problems  discussed.  In  addition  a  brief  overview 
of  the  more  advanced  wavelet  modulus  maximum  method  (WTMM)  presented.  A  segmentation  method 
is  present  which  further  increases  the  accuracy  of  the  WTMM  method. 


2.1  Box-counting 

Traditionally,  an  objects  fractal  dimension  can  be  calculated  by  covering  it  with  boxes  of  varying  size  e.  The 
fractal  dimension  is  then  given  by  the  scaling  relation: 

N (e)  oc  e~D ,  (1) 

where  D  is  the  fractal  dimension,  and  N(e)  is  the  number  of  boxes  of  size  e.  From  this,  it  is  straight  forward 
to  show  that  a  line  has  a  fractal  dimension  of  one  as  the  number  of  boxes  required  scales  linearly,  similarly 
for  a  square  the  scaling  relation  is  quadratic  and  has  a  fractal  dimension  of  two. 

Multifractals,  like  fractals,  relate  the  number  of  boxes  N(e)  to  the  size  of  each  box  e  and  are  characterised 
by  a  similar  equation: 

N(e)  oce-/(a),  (2) 

where  instead  of  a  single  exponent,  there  is  a  spectrum  of  exponents  /(a),  each  with  a  relative  strength  or 
significance  a. 

An  overview  of  the  box-counting  method  as  applied  to  solar  line-of-sight  magnetogram  data  is  provided 
in  (Conlon  et  al.,  2008).  Briefly,  the  main  steps  are  to  calculate  the: 

1.  magnetic  field  with  each  box  i  of  size  e,  Be(i). 

2.  partition  function  for  each  box  size  e,  ze{q). 

3.  normalised  flux  for  eacbox  i  of  size  e,  Be(i). 
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Fig.  1 :  Left:  MDI  magnetogram  taken  28  October  2003.  Center:  Wavelet  modulus  used  in  the 
fractal  calculations  with  no  segmentation.  Right:  Wavelet  modulus  used  in  the  fractal  cal¬ 
culations  following  segmentation  Kestener  et  al.  (2009). 

4.  strength  or  significance  of  a  measure,  a(q). 

5.  distribution  of  points  of  a  given  concentration  strength,  /(a). 

/(a)  for  q  =  0  corresponds  to  the  Hausdroff  or  fractal  dimension. 

Due  to  the  discrete  nature  of  the  box-counting  method,  it  has  many  associated  errors  (Conlon  et  al.,  2008; 
Alber  &  Peinke,  1998).  Methods  have  been  developed  to  overcome  these  discritisation  errors  and  results 
show  they  can  successfully  return  the  entire  or  part  of  the  multifractal  spectrum.  However,  these  methods 
tend  to  be  finally  tuned  for  the  specific  data  sets  used. 

In  addition,  there  is  a  added  error  due  to  possible  mixing  of  fractal  and  multifractal  spectrums,  either  due 
to  noisy  data  or  the  embedding  of  the  system  in  question  within  another  surrounded  system.  Previously 
this  error  was  reduced  by  limiting  the  field  of  view  to  only  cover  the  object  in  question.  However  given  the 
fractal  geometry  of  the  object  in  question  these  results  can  vary. 


2.2  Wavelet  Modulus  Maximum  Method  (WTMMM) 

In  order  to  overcome  the  discrete  nature  of  the  box-counting  method  a  wavelet  based  method  was  devel¬ 
oped.  The  wavelet  modulus  maximum  method  (WTMM)  method  replaces  the  boxes  from  the  traditional 
box-counting  method  with  wavelets.  Wavelets  have  the  advantage  of  been  defined  for  both  finite  and  dis¬ 
crete  domains,  and  in  term  as  a  fuzzy  boxes.  The  wavelet  based  method  however  returns  the  so-called 
singularity  spectrum,  D(h).  This  differs  from  the  multifractal  spectrum,  /(a)  as  follows: 

f(a)  =  D(h),  (3) 

a  =  h  +  Dimension(image),  (4) 

where  h  is  the  Holder  exponent.  In  order  to  calculate  the  D{h)  singularity  spectrum  a  2D  WTMM  method 
developed  by  Arneodo,  A.,  &  Decoster,  N.  (2000)  is  used.  An  overview  of  the  WTMM  method  is  provided 
in  Kestener  et  al.  (2009).  Briefly  the  main  steps  of  the  WTMM  method  are: 

1.  Compute  of  the  2D  continuous  wavelet  transform  of  the  input  image. 

2.  Extract  the  WTMM  edges  at  each  scale. 

3.  Calculate  the  local  maximum  along  the  wavelet  transform  (WT)  skeleton  and  their  linking  across 
each  scale. 

4.  Calculate  the  partition  functions  from  the  WT  skeleton. 

5.  Calculate  the  remaining  fractal  and  multifractal  parameters. 
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As  mentioned  already,  another  major  source  of  error  in  the  calculation  of  fractal  and  multifractal  properties 
is  the  contribution  from  alien  process  altering  the  true  spectrum.  In  order  to  overcome  this,  a  segmentation 
in  wavelet  modulus  space  is  proposed.  This  methods  requires  that  the  two  process,  in  this  case  active 
and  quite  Sun  magnetic  features,  are  significantly  separate  and  distinct  from  each  other  in  some  wavelet 
modulus  space.  Kestener  et  al.  (2009)  found  this  for  a  quite  and  active  region  magnetic  field  and  proposed 
at  cut  off  at  a  modulus  strength  of  40.  The  segmentation  method  is  shown  to  greatly  increase  the  accuracy 
of  the  returned  multifractal  spectrum  for  both  test  and  real  data.  An  example  of  the  segmentation  methods 
ability  to  selectively  extract  wavelet  modulus  related  to  active  region  is  shown  in  Figure  1 .  As  is  shown  the 
segmentation  method  selects  wavelet  modulus  related  to  the  active  regions  feature  on  the  solar  disk  and 
removes  weaker  modulus  related  to  the  quiet  Sun. 


3  Results 


This  work  represent  an  overview  of  a  much  larger  body  of  work  (Conlon  et  al.,  2009;  Kestener  et  al.,  2009), 
which  has  shown  the  WTMM  and  the  associated  segmentation  method  to  be  a  stable  and  accurate  math¬ 
ematical  tool  for  the  analysis  of  magnetic  fields  on  the  solar  disk.  Two  examples  are  now  described  in 
detail.  These  example  further  highlight  the  proposed  thresholds  of,  D{h)  >  1.2  and  h  >  —0.7  as  indicator 
of  increased  flaring  potential. 


3.1  NOAM  0488 

NOAA  10488  emerged  onto  the  solar  disk  on  the  26  October  2003  and  continued  to  grow  rapidly  until  the 
29  October  2003  (Figure  2).  The  ability  of  the  segmentation  methods  to  detect  the  presense  of  active  regions 
in  quiet  Sun  data  is  highlight  by  the  rise  in  D{h)  and  h  prior  to  the  increase  in  the  other  physical  parameters. 
Prior  to  the  regions  emergence  only  a  limited  number  of  wavelet  modulus  moments  exceed  the  threshold, 
resulting  in  D(h )  and  h  values  that  would  characterise  poorly  resolved  dusty  data.  As  the  regions  forms 
and  develops  into  a  coherent  structure  the  fractal  dimension  and  Holder  exponent  rise  steadily. 

While  on  disk  NOAA  10488  underwent  two  flaring  periods.  The  first  occurred  over  an  eighteen  hour  period 
starting  at  14:02  UT  on  the  27  October  2003  and  contained  five  C-class  flares  and  a  Ml. 9  flare.  The  second 
started  at  18.10  UT  on  the  29  October  2003,  lasted  for  two  days  and  resulted  in  8  C-class  flares  and  two  M- 
class  flares.  During  the  regions  evolution  the  associated  fractal  dimension,  D(h),  is  seen  to  increase  above 
the  proposed  threshold  for  flaring  (McAteer,  Gallagher,  and  Ireland,  2005).  In  addition  flaring  periods  are 
also  associated  a  Holder  exponent,  h,  greater  than  the  proposed  threshold  of  Conlon  et  al.  (2009).  Given 
that  the  Holder  exponent  relates  the  continuity  of  the  underlying  magnetic  field,  a  Holder  exponent  close 
to  0  represents  magnetic  features  characteristic  of  large  gradient  or  step-like-functions. 


3.2  NOAA  10953 

NOAA  10953  rotated  onto  the  solar  disk  on  the  25  April  2007  already  fully  developed.  While  on  disk  NOAA 
10953  underwent  no  flaring  periods  during  the  reported  observation.  Figure  3  highlights  the  changing 
fractal  dimension  and  Holder  exponent  of  NOAA  10953.  At  no  time  during  the  observation  period  was 
both  the  regions  fractal  dimension  and  Holder  exponent  greater  than  the  proposed  threshold  from  Conlon 
et  al.  (2009).  Compared  with  the  results  from  Conlon  et  al.  (2009)  this  is  further  evidence  that  the  proposed 
threshold  are  a  valid  indicator  for  a  regions  ability  to  flare. 


4  Future  Direction 


The  results  represented  here  and  in  Conlon  et  al.  (2009);  Kestener  et  al.  (2009)  are  related  to  a  single  point 
on  the  multifractal  spectrum,  q  =  0.  The  study  of  the  remaining  parts  of  the  multifractal  spectrum  is  left  for 
a  later  date. 


4  Future  Direction 
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These  results  further  highlight  the  ability  of  multiscale  methods  to  detect  characteristics  changes  in  the 
make-up  of  elements  with  active  region  on  the  photosphere.  Together  with  other  physical  parameters,  the 
multifractal  spectrum  may  make  it  possible  to  detect  conditions  favourable  to  flaring  in  an  active  region. 
Further  work  is  required  to  validate  and  constrain  any  thresholds  that  may  or  may  not  exist  in  reality. 


Area  (Mm2) 


4  Future  Direction 
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25-Oct-2003  19:11  27-Oct-2003  07:59  28-Oct-2003  11:11  29-Oct-2003  12:51 


26-0ct  27-0ct  28-0ct  29-0ct  30-0ct  31-0ct 


Fig.  2:  The  evolution  of  NOAA  10488.  (Top  panel)  MDI  magnetogram  image  for  NOAA  10488  on 
25  October  2003  at  19:11  UT,  27  October  2003  at  07:59  UT,  28  October  2003  at  11:11  UT  and  29 
October  2003  at  12:51  UT.  (Second  panel)  D{h)  for  q  =  0.  (Third  panel)  h  for  q  =  0.  (Fourth 
panel)  The  total  unsigned  magnetic  flux  (Mx).  (Bottom  panel)  Total  area  (Mm2).  Associated 
C-class  flares  are  indicated  by  thin  arrows;  bolder  arrows  indicate  M-class  flares.  Vertical 
dotted  lines  indicate  the  time  of  MDI  magnetograms  in  the  regions  evolution.  Dashed  hori¬ 
zontal  line  in  the  first  panel  highlight  a  previously  proposed  threshold  for  flaring  (McAteer, 
Gallagher,  and  Ireland,  2005).  Dashed  horizontal  line  in  the  second  panel  shows  our  pro¬ 
posed  threshold.  Conlon  et  al.  (2009) 


Area  (Mmz) 
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Fig.  3:  The  evolution  of  NOAA  10953.  (Top  panel)  MDI  magnetogram  image  for  NOAA  10953  on 
30  April  2007  at  12:51  UT,  1  May  2007  at  14:27  UT,  and  2  May  2007  at  04:51  UT.  (Second 
panel)  D(h )  for  q  =  0.  (Third  panel)  h  for  q  =  0.  (Fourth  panel)  The  total  unsigned  mag¬ 
netic  flux  (Mx).  (Bottom  panel)  Total  area  (Mm2).  Dashed  horizontal  line  in  the  first  panel 
highlight  a  previously  proposed  threshold  for  flaring  (McAteer,  Gallagher,  and  Ireland, 
2005).  Dashed  horizontal  line  in  the  second  panel  shows  our  proposed  threshold. 


4  Future  Direction 


8 


References 

Abramenko,  V.  L:  2005,  Astro.  Phys.  J.,  629, 1141 

Abramenko,  V.  I.:  2005,  Solar  Phy.,  228, 29 

Abramenko,  V.,  Yurchyshyn,  V.,  &  Wang,  H.  2008,  ApJ,  681, 1669 

Alber,  ML,  &  Peinke,  J.  1998,  Phys.  Rev.  E,  57, 5489 

Arneodo,  A.,  Decoster,  N.,  &  Roux,  S.  G.  2000,  Eur.  Phys.  J.  B,  15,  567 

Aschwanden,  Ml.  J.,  &  Aschwanden,  P.  D.  2008,  ApJ,  674, 530 

Conlon,  P.  A.,  Gallagher,  P.  T.,  McAteer,  R.  T.  J.,  Ireland,  J.,  Young,  C.  A.,  Kestener,  P,  Hewett,  R.  J.,  & 
Maguire,  K.  2008,  Solar  Phys. ,  248, 297 

PA.  Conlon  ,  P.  Kestener  ,  R.T.J.  McAteer  ,  PT.  Gallagher  ,  L.  Fennell ,  A.  Khalil ,  and  A.  Arneodo,  2009, 
A&A,  in  prep. 

Gallagher,  P.  T.  et  al:  2006,  Astrophysics  And  Space  Science  Library,  344, 1.2 

Georgoulis,  M.  K.:  2005,  Solar  Phy.,  228, 5 

Georgoulis,  M.  K.  2008,  Geophysical  Research  Letters,  35,  6 

Hewett,  R.  J.,  Gallagher,  P.  T.,  McAteer,  R.  T.  J.,  Young,  C.  A.,  Ireland,  J.,  Conlon,  P.  A.,  &  Maguire,  K.  2008, 
Solar  Phys.,  248, 311 

Kestener,  R,  et  al.  2009,  A&A,  in  prep. 

Mandelbrot,  B.  B.  1974,  Journal  of  Fluid  Mechanics,  62, 331 

Mandelbrot,  B.  B.  1977,  San  Francisco:  W.H.  Freeman,  1982,  Revised  edition  of:  Fractals  (1977), 

McAteer,  R.  T.  J.,  Gallagher,  P.  T.,  &  Ireland,  J.  2005,  ApJ,  631,  628 
Miesch,  M.  S.  2005,  Living  Reviews  in  Solar  Physics,  2, 1 
Schroeder,  M.  1991,  New  York:  Freeman,  1991, 


5  Conference  presentation/Publication 


9 


5  Conference  presentation/Publication 


Reference  Papers:  Turbulence,  Complexity,  and  Solar  Flares,  R.T.J.  McAteer,  P.T.  Gallagher,  P.A.  Conlon,  Ad¬ 
vances  in  Space  Science  Research,  submitted. 

Characterising  complexity  in  Compound  Systems  II:  Magnetic  Fields,  Flares  &  Forecasts.  P.A.  Conlon ,  P.  Kestener 
,  R.T.J.  McAteer  ,  P.T.  Gallagher  ,  L.  Fennell ,  A.  Khalil ,  and  A.  Arneodo,  Astronomy  &  Astrophysics,  in 
prep. 

Characterising  complexity  in  Compound  Systems  I:  Segmentation  in  Wavelet-Space.  P.  Kestener  ,  P.A.  Conlon ,  A. 
Khalil ,  A.  Arneodo ,  R.T.J.  McAteer ,  P.T.  Gallagher ,  and  L.  Fennell,  Astronomy  &  Astrophysics,  in  prep. 

Multifractal  Properties  of  Evolving  Active  Regions  P.  A.  Conlon,  P.  T.  Gallagher,  R.  T.  J.  McAteer,  J.  Ireland,  C. 
A.  Young,  P.  Kestener,  R.  Heweitt,  K.  Maguire.  Solar  Phys.,  2, 4.8, 297, 2008. 

Book  Chapter: 

Solar  Activity  Monitoring  Gallagher,  P.  T.,  McAteer,  R.  T.  J.,  Young,  C.  A.,  Ireland,  J.,  Hewett,  R.,  Conlon,  P, 
In  J.  Lilensten  (ed)  "In  SPACE  WEATHER:  Research  towards  Applications  in  Europe",  Springer,  (2006). 

Poster  paper: 

Multifractal  analysis  of  magnetic  fields  in  the  solar  atmosphere  Conlon,  P,  Gallagher,  P.  T.,  McAteer,  R.  T.,  J., 
Young,  C.  A.,  Ireland,  J.,  K.  Maguire  European  Space  Weather  Week  III.  Brussels,  Belgium.  November  2006. 

Multifractal  Analysis  of  Sunspot  Magnetic  Fields  Conlon,  P.,  Gallagher,  P.  T.,  McAteer,  R.  T.,  J.,  Young,  C.  A.,  Ire¬ 
land,  J.,  K.  Maguire  RAS  Specialist  Discussion  Meeting,  Solar-Terrestrial  Magnetic  Fields  &  Reconnection. 
London,  England  September  2006. 

Multifractal  analysis  of  magnetic  fields  in  the  solar  atmosphere  Conlon,  P,  Gallagher,  P.  T.,  McAteer,  R.  T.,  J., 
Young,  C.  A.,  Ireland,  J.,  American  Astronomical  Society,  Durham,  New  Hampshire,  USA.  June  2006. 

Talks: 

Magnetic  Complexity  in  Solar  Active  Regions  R.T.J.  McAteer,  Invited  review  NASA  All  Clear  Workshop,  Boul¬ 
der,  CO,  USA,  April  2009. 

Turbulence,  Complexity,  and  Solar  Flares,  R.T.J.  McAteer,  P.T.  Gallagher,  P.A.  Conlon,  European  Space  Weather 
Week,  Brussels,  Belgium,  November  2008. 

Multifractal  analysis  of  solar  magnetic  fields:  Implications  for  space  weather  forecasting  Conlon,  P.  Contributed 
talk  at  Solar  Image  Processing  Workshop  III,  Trinity  College  Dublin,  Dublin,  Ireland,  September  2006. 


